MSA 

TP 

1120 

o.l 

NASA Technical Paper 1120 


Loan copy: return -n S » 

AFWL TECHNICAL LIftRA f = 

KIRTLAND AFB, N'. m J ^ = 

Numerical Study of Transonic Flow j 
Over Oscillating Airfoils Using 
the Full Potential Equation 


■ L \ 

Koji Isogai 


APRIL 1978 


\ •• S t 1 J 7 / -v. 
sV- • \ 

’V.F" "\ J \ 

/ 

5 Y v 

' i 'i ’;1|- \ 

-j * . *9?S ■: 


4# £?*•■*«*, 


:\ 

:'\ -cv 




TECH LIBRARY KAFB, NM 


TECH LIBRARY KAFB, NM 

■iiiiini 

□lBMMTD 

NASA Technical Paper 1120 


Numerical Study of Transonic Flow 
Over Oscillating Airfoils Using 
the Full Potential Equation 


Koji Isogai 

Langley Research Center 
Hampton, Virginia 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
information Office 


1978 



SUMMARY 


The behavior of unsteady aerodynamic loadings on airfoils oscillating in 
transonic flow has been investigated numerically with particular attention given 
to supercritical airfoil sections. A previously developed finite difference 
method, which is based on the full potential equation and which uses a quasi- 
conservative scheme for proper capture of shock wave motion, was employed for 
the present study. The unsteady aerodynamic pressure and load distributions on 
several different airfoil sections are presented with particular emphasis on 
the effects of free-stream Mach number, reduced frequency, and mean angle of 
attack. These parameters are demonstrated to have a significant effect on the 
behavior of the unsteady aerodynamic loadings. Comparisons of the present cal- 
culations with the exact inviscid solution and with the experimental results 
are also presented. 


INTRODUCTION 

Recent development of supercritical airfoil sections has led to improve- 
ments in aerodynamic efficiency at high subsonic speeds. This situation has 
stimulated considerable interest in unsteady aerodynamics of supercritical air- 
foil sections (refs. 1 and 2). This interest was reinforced when wind-tunnel 
tests (ref. 3) showed that a supercritical wing experiences a more pronounced 
dip of the flutter boundary in the transonic regime in comparison to a struc- 
turally similar conventional wing. 

Recently, reference 4 presented a method for solving the full potential 
equation by a time marching finite difference technique. The method was a quasi- 
conservative scheme to capture shock wave motion properly. The capability of 
the method for treating a supercritical airfoil section was well demonstrated 
by numerical examples in reference 4. No systematic investigation of the behav- 
ior of unsteady aerodynamic loadings was attempted, however. Since the super- 
critical transonic flow field (especially about supercritical airfoil sections) 
is sensitive to changes of free-stream Mach number and angle of attack, a con- 
siderable effect of these parameters on the behavior of unsteady aerodynamic 
loadings is expected. In this paper, the method of reference 4 is applied in 
an investigation of the behavior of unsteady aerodynamic loadings on oscillating 
airfoils (including supercritical airfoils). Particular emphasis is given to 
the effects of mean angle of attack, free-stream Mach number, reduced frequency, 
and airfoil shape. 

The items investigated in the present paper can be summarized as follows: 

(a) Behavior of unsteady aerodynamic loadings on supercritical airfoils at 
design and off-design conditions. 



(b) Comparison of unsteady aerodynamic forces, pitching moments, and load 
distributions between a supercritical airfoil and a conventional airfoil. 

(c) Effect of mean angle of attack on unsteady load distributions on a 
supercritical airfoil. 

(d) Behavior of higher harmonics of unsteady aerodynamic loadings. 


SYMBOLS 


3 n 


harmonic component of total lift coefficient or pitching-moment coef- 
ficient (n = 0, 1, . . . , N) 




harmonic component of pressure loading coefficient (n = 0, 1, 


, N) 


semi chord 
lift coefficient 

first harmonic of lift coefficient per unit 0 

pitching-moment coefficient taken about midchord (positive nose up) 
first harmonic of pitching-moment coefficient per unit 0 


pressure coefficient, 


P - Poo 

1 P_V 2 


ACn 



pressure coefficient at sonic condition 
pressure loading coefficient ^ACp = Cp - Cp) 


in-phase component of the first harmonics of the pressure loading 
coefficient per unit 0 and 6, respectively 


out-of-phase component of the first harmonics of the pressure loading 
coefficient per unit 0 and 6, respectively 


= vPT 
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reduced frequency, boo/V^ 
lift 

Mach number 

highest harmonic 

nth harmonic (1, 2, 3, . . . , N) 

pressure 

time, seconds 

dimensionless time, (V^/bjT 
total velocity 

coordinate in chordwise direction normalized by b and measured from 
midchord 

angle of attack 

mean angle of attack 

angle of control surface (flap) deflection 
amplitude of pitching oscillation, radians 
density 

phase angle between total lift and pitching moment, respectively, and 
displacement in pitching oscillation, positive for total lift or 
pitching-moment leading displacement 


circular frequency 
j | magnitude 

Subscripts: 

t stagnation condition 

00 free-stream condition 

Superscripts : 

U,L upper and lower surface of airfoil 
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RESULTS FOR A NACA 64A006 AIRFOIL WITH AN OSCILLATING FLAP 


Comparisons between the present full potential calculation and the invis- 
cid results obtained by using the Euler equations in reference 5 (Magnus and 
Yoshihara) are presented in this section. The case considered is a NACA 64A006 
airfoil with oscillating quarter-chord flap at M^ = 0.875 and k □ 0.234. 
Comparisons of the steady pressure distributions for the two calculations 
are shown in figure 1. For reference, experimental results from reference 6 
(Tij deman and Schippers) are also shown. The agreement between the present full 
potential calculation and the exact inviscid result is satisfactory. The dis- 
crepancy between experiment and calculations can possibly be attributed to the 
wind-tunnel wall interference and boundary-layer effects as discussed in 
reference 7. 

Comparisons of unsteady pressure distributions between the two calculations 
are shown in figures 2(a) to 2(d) for upper and lower surface pressures at two 
different phases of the oscillation of the quarter-chord flap. Although the 
present full potential calculation predicts a slightly stronger shock wave than 
that of the exact result, the overall agreement between the two calculations is 
good. The prediction of the stronger shock wave by the present method is rea- 
sonable because the shock jump relation of the full potential equation gives a 
stronger shock wave than that of the Rankine-Hugoniot shock jump relation 
(ref. 8). The same degree of agreement between the two calculations was 
obtained for other phases of the flap oscillation. 

In-phase and out-of-phase components of the unsteady load distributions 
for the same case are compared with experimental results obtained from refer- 
ence 6 in figures 3(a) and 3(b). Although qualitative agreement between the 
calculation and experiment is good, there is some quantitative discrepancy. 
Reference 7 pointed out that the discrepancy between calculations, including 
the boundary-layer correction and experimental results, can be greatly reduced 
if the nominal free-stream Mach number of the experiment (0.875) is reduced to 
0.850 to allow for the wind-tunnel wall interference. 

Although the present computer code is not capable of treating the boundary- 
layer effect, the discrepancy between the calculation and experiment in fig- 
ures 3(a) and 3(b) can be reduced by adjusting the free-stream Mach number (of 
the calculation) so that the theoretical shock position of the mean steady flow 
comes as close as possible to the experimental shock position obtained at 
M^ = 0.875. The closest agreement of the steady shock position between the 
theory and experiment is obtained at M^ = 0.860 rather than at M^ = 0.850 
obtained in reference 7, as shown in figure 4. The corresponding unsteady 
results are compared with the experimental results from reference 6 in fig- 
ures 5(a) and 5(b). The agreement between the two results is markedly improved. 
This result indicates that the mismatch of the shock position between the calcu- 
lation and experiment is the predominant factor in the disagreement of the load 
distribution between the theory and experiment shown in figures 3(a) and 3(b). 
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RESULTS FOR A SUPERCRITICAL AIRFOIL OSCILLATING IN PITCH 


AT DESIGN AND OFF-DESIGN CONDITIONS 

Calculations of steady pressure distributions and unsteady load distribu- 
tions were performed for the 70-10-13 supercritical airfoil of reference 9 for 
three different, flow conditions: 

Case 1: = 0.71, a m = 0° (off design) 

Case 2: = 0.70, (% = 0° (design) 

Case 3: = 0.68, a m = 1° (off design) 

The steady pressure distribution for each case is shown in figure 6(a). The 
calculations for unsteady load distributions corresponding to these steady flow 
conditions were performed for the airfoil oscillating in pitch about the mid- 
chord axis with an amplitude of 1° at the reduced frequency of 0.10. The 
unsteady load distribution (the first harmonic) for each case is shown in fig- 
ure 6(b). The in-phase and out-of-phase components are shown by the solid line 
and the dashed line, respectively. 

Examination of variation of the instantaneous pressure distributions^ dur- 
ing the cycle of oscillation for these three cases reveals that the patterns of 
the unsteady shock wave motion on the upper surface are very different for each 
case. In case 1, the relatively strong shock wave oscillates around its mean 
position with relatively small excursion amplitude (about 4 percent of the 
chord) during the cycle of the oscillation. In case 2, a relatively weak shock 
wave appears and disappears during the cycle of the oscillation in the neighbor- 
hood of the part of the airfoil where the local supersonic region of the mean 
steady flow terminates. In case 3, as discussed later in detail, the relatively 
weak shock wave on the upper surface oscillates back and forth with large excur- 
sion amplitude over the forward part of the airfoil during the cycle of the 
oscillation. 

In figure 6(b), the differences in the unsteady shock wave motion between 
these three cases are clearly reflected in the unsteady load distributions. In 
case 1 , the relatively strong shock wave causes a peak value near the three- 
quarter-chord point for both the in-phase and out-of-phase components whereas 
such a peak value caused by the shock wave is not observed in the load distri- 
bution for case 2. In case 3, a bulge appears in the pressure distribution on 
the front part of the airfoil. Comparison of the full potential unsteady load 
distribution for case 3 with that of flat-plate theory (see fig. 7) indicates 
that this bulge contributes considerably to the unsteady total lift and pitching- 
moment derivatives. Examination of the variation of instantaneous pressure dis- 
tributions during a cycle of oscillation, shown in figures 8(a) to 8(d), reveals 
that the formation of this bulge is closely associated with a shock wave on the 
upper surface which oscillates back and forth with relatively large excursion 
amplitude. The formation of such a bulge in the unsteady (upper surface) pres- 

^The instantaneous pressure distributions are shown in this paper only for 
case 3 (see fig. 8). 
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sure distribution was also observed experimentally in reference 2 on a NLR 7301 
"shock free” airfoil oscillating in pitch at the design condition. 


COMPARISON OF RESULTS FOR SUPERCRITICAL AND CONVENTIONAL 
AIRFOILS OSCILLATING IN PITCH 

According to the recent experimental study of reference 3, a wing with a 
supercritical airfoil experienced a more pronounced dip of a flutter boundary 
than a structurally similar conventional wing. Unsteady aerodynamic forces are 
clearly responsible for this phenomenon. Therefore, the difference of the 
unsteady aerodynamic derivatives and aerodynamic loadings between the two air- 
foils was investigated. For this purpose a 79-03-12 supercritical airfoil 
(ref. 9, Bauer et al.) and a NACA 0012 conventional airfoil were selected. 

Both airfoils have the same ratio of thickness to chord. In figure 9, the 
steady pressure distribution calculated on the 79-03-12 airfoil for the design 
condition of = 0.79 and a = 0° is compared with the design pressure dis- 
tribution from reference 9. The agreement between the present calculation and 
the design pressure distribution is satisfactory. In figure 10, the steady 
pressure distributions of each airfoil at three different free-stream Mach num- 
bers for an angle of attack of 0° are shown. The unsteady calculations have 
also been performed for these three Mach numbers. The airfoils were considered 
to be oscillating in pitch about the midchord axis with an amplitude of 1° at a 
reduced frequency of 0.10. 

The magnitudes and phase angles of the unsteady lift (first harmonic) and 
pitching moment (first harmonic) for both airfoils are shown in figures 11(a) 
and 11(b), respectively. With the exception of the M^ □ 0.75 case, the mean 
angle of attack was 0° for all calculations. At M^ = 0.75, results were also 
computed for an angle of attack of 1° to illustrate the effect of the mean angle 
of attack. The values predicted by the flat-plate theory are also presented in 
the same figures for comparison. The unsteady load distributions (in-phase and 
out-of-phase components of the first harmonic), which correspond with the steady 
flows shown in figure 10, are compared in figures 12(a) to 12(c) at each Mach 
number. For the unsteady lift (see fig. 11(a)), the difference between the two 
airfoils is surprisingly small. However, this apparent agreement of the total 
lifts is rather fortuitous. As shown in figures 12(a) to 12(c), there is a con- 
siderable difference in the load distribution between the two airfoils except 
those of M^ = 0.70 for which the flows are almost subcritical. This dissimi- 
larity in the load distribution between the two airfoils is reflected by the 
difference in pitching moment (about midchord axis), which becomes greater as 
the Mach number is increased as shown in figure 11(b). It is also observed in 
figures 11(a) and 11(b) that the change of the mean angle of attack affects the 
unsteady lift and moment differently between the two airfoils. In figure 12(c) 
there is still a considerable nonlinear thickness effect on total lift and 
pitching moment at M^ = 0.70 even though the flow is almost subcritical. 

This effect can be confirmed by comparing load distributions for the nonlinear 
results with flat-plate theory. 
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RESULTS FOR EFFECT OF MEAN ANGLE OF ATTACK FOR A 


SUPERCRITICAL AIRFOIL OSCILLATING IN PITCH 

Mean angle of attack was shown to affect both the unsteady total lift and 
pitching moment for the 79-03-12 and NACA 0012 airfoils (figs. 11(a) and 11(b)). 
Angle-of-attack effects are, therefore, examined more closely for the 79-03-12 
airfoil. In figure 13, steady pressure distributions for the 79-03-12 super- 
critical airfoil at M^ = 0.75 at two different angles of attack (0° and 1°) 
are shown. At a = 0°, the steady flow is slightly supercritical near the 
leading-edge region without shock wave, whereas at a = 1° the local supersonic 
region on the upper surface is increased, being terminated with a weak shock 
wave near the 40 -percent-chord point. 

For unsteady calculations the airfoil is oscillated in pitch with an ampli- 
tude of 1° at k = 0.10 about these different mean angles of attack (0° and 1°). 
The in-phase and out-of-phase components (the first harmonic) of the unsteady 
load distributions corresponding to these two different mean angles of attack 
are compared in figure 14. The mean angle of attack has a striking effect. 

Such a difference in the load distribution comes from the difference of the 
unsteady shock wave pattern in the two cases. For the case of ot m = 0°, the 
shock wave does not appear during the cycle of oscillation, whereas for the 
case of a = 1°, the relatively weak shock wave oscillates back and forth over 
the forward portion of the upper surface of the airfoil. This behavior is evi- 
dent in the variation of the instantaneous pressure distribution shown in 
figure 15. 


RESULTS FOR HIGHER HARMONICS OF LOADING 
ON AIRFOILS OSCILLATING IN PITCH 

In transonic flow, effects of the shock wave displacement cause nonsinusoi- 
dal variations in the local pressures on some parts of the airfoil. In this 
section, the chordwise distributions of higher harmonic components of unsteady 
aerodynamic loadings on two different airfoils which exhibit different unsteady 
shock wave pattern, are examined. The first case considered is the 79-03-12 
supercritical airfoil oscillating in pitch about the midchord axis at M^ = 0.75 
and k = 0.10. The mean angle of attack and the amplitude of the oscillation 
are both 1°. The steady pressure distribution for this case is shown in fig- 
ure 13(b). The instantaneous pressures for the oscillating case are shown in 
figures 15(a) to 15(d). The relatively weak shock wave oscillates around the 
40-percent-chord point with relatively large shock excursion amplitude (about 
25-percent chord). The shock wave changes its strength considerably during the 
cycle of oscillation (the shock wave disappears at its maximum upstream position 
as seen in fig. 15(d)). 

Harmonic analysis of the time-varying load distributions was made. The 
resulting coefficients of the harmonics are used to define AC p in the follow- 
ing equation: 
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N 

ACp = ag + 6 £a n cos (nkt) + b n sin (nkt)J 

n=1 

The chordwise distributions of b n (n = 1, 2, 3) and a n (n = 1, 2, 3) are 
plotted in figures 16(a) and 1 6 ( b ) , respectively. In these figures, two char- 
acteristics of the higher harmonics (n = 2, 3) are evident. First, the region 
where the magnitudes of the higher harmonics are appreciable is confined to the 
portion of the airfoil over which the shock wave oscillates back and forth (see 
figs. 15(a) to 15(d)). Second, the higher harmonics exhibit a wavy pattern in 
the chordwise direction. As the simplified model used in reference 2 shows, 
the existence of the wavy pattern of the higher harmonics is the natural conse- 
quence of the excursion of the sharp pressure rise associated with the moving 
shock wave. 

To determine the effect of reduced frequency on harmonics of the load dis- 
tributions, the unsteady calculation for k = 0.30 was also performed without 
changing the other flow conditions. The results are shown in figures 17(a) and 
17(b). Some effect of the reduced frequency is observed on both the first har- 
monic and higher harmonics. The region where the higher harmonics are relatively 
large is decreased for the higher reduced frequency. This result is caused by 
a reduction in shock excursion amplitude at the higher reduced frequency. 

The second case considered is that of a NACA 0012 airfoil oscillating in 
pitch about midchord axis at = 0.79 and k = 0.10. The steady pressure 

distribution and the variation of the instantaneous pressure distribution are 
shown in figure 18 and figures 19(a) to 19(d), respectively. The relatively 
strong shock wave oscillates about the midchord position with a relatively small 
excursion amplitude (about 12-percent chord). The chordwise distributions of 
harmonic components b n (n = 1, 2, 3) and a n (n = 1, 2, 3) of the loading are 
shown in figures 20(a) and 20(b), respectively. Again, the region where the 
higher harmonics are appreciable is confined to the part of the airfoil over 
which the shock wave excursion occurs. The higher harmonics differ from the 
previous case because the magnitude of the second and third harmonics (b2, bj, 
a2, and 83) is relatively" small in comparison with the very large peak values 
of the first harmonics (b-| and a-j) observed in the region near the oscillat- 
ing shock wave. 

As pointed out in reference 2, the contribution of higher harmonic compo- 
nents to the unsteady total lift and pitching moment might be relatively small 
in comparison with the contribution of the first harmonic because of their wavy 
chordwise distributions. This result can be confirmed in table I where harmonic 
components of the total lift and pitching moment for both cases (79-03-12 
supercritical airfoil and NACA 0012) are shown. The harmonic components A n 
(n = 0, 1,2, 3) and B n (n = 1, 2, 3) of the total lift and pitching moment 
are defined by the following equation: 

C L N 

or = Aq + 0 [A n cos (nkt) + B n sin (nkt)J 

^m n=1 
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The magnitude of the higher harmonic components is less than 8 percent of the 
first harmonic components for case 1 and less than 2 percent of the first har- 
monic components for case 2. 

CONCLUDING REMARKS 

A numerical study has been performed to investigate the behavior of the 
unsteady aerodynamic forces and moments, load distributions, and pressure dis- 
tributions on airfoils (including supercritical airfoils) oscillating in a 
transonic flow. Particular emphasis has been placed on the effect of free- 
stream Mach number, mean angle of attack, and reduced frequency. It has been 
shown that these parameters have a significant effect on the behavior of the 
unsteady aerodynamic loadings. Generally, the change of the parameters pro- 
duced a different pattern of the unsteady shock wave motion, which played the 
predominant role in determining the unsteady aerodynamic load distributions. 

Since the present results are based on the inviscid flow calculations, the 
potentially important effect of the boundary layer has been neglected. However, 
good agreement for the unsteady load distributions between calculation and exper 
iment was obtained by adjusting the free-stream Mach number so that the theoreti 
cal and experimental steady shock positions are in agreement. Thus, the invis- 
cid flow appears to play a dominant role in determining the behavior of the 
unsteady aerodynamic loadings as long as the boundary layer remains attached. 
Although the inclusion of the boundary-layer effect is obviously necessary, the 
present results, based on the inviscid flow calculations, give insight into the 
unsteady transonic aerodynamics of oscillating airfoils. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 17, 1978 
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TABLE I.- HARMONIC COMPONENTS OF UNSTEADY TOTAL LIFT AND PITCHING MOMENT 
ON OSCILLATING 79-03-12 SUPERCRITICAL AIRFOIL AND NACA 0012 AIRFOIL 


(a) 79-03-12 supercritical airfoil; M = 0.75; k = 0.10; 

a m = 1°; 0 = 1° 



A 0 

Al 

Bi 

a 2 

b 2 

a 3 

b 3 

c L 

0.4872 

-4.0362 

5.5908 

0.2429 

0.0104 

-0.0331 

0.0135 

Cm 

.0677 

-1.4237 

1.1547 

.0645 

.0915 

.0094 

.0149 


(b) 

NACA 0012 airfoil 

I M co = 1 

3.79; k 

= 0.10; 

<% = o°; 

o 

t — 

II 

CD 


A 0 

A 1 

Bi 

a 2 

b 2 

a 3 

b 3 

Cl 

-0.0019 

-3.4010 

6.5862 

0.0066 

-0.0399 

-0.0099 

-0.0262 

Cm 

-.0002 

-.9547 

1.1693 

.0028 

-.0202 

.0095 

-.0123 
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Magnus and Yoshihara (ref. 5) 

(Euler equations) 

Present calculation (full potential) 



Figure 1.- Steady pressure distributions for NACA 64A006 airfoil. 

M = 0.875; a = 0°. 

oo ' 




(a) Upper surface; kt = 30°. 



(b) Lower surface; kt = 30°. 


Figure 2.- Unsteady pressure distributions on NACA 64A006 airfoil 
with oscillating quarter-chord flap. M = 0.875; k = 0.234. 





(a) In-phase component. 



(b) Out-of-phase component. 

Figure 3.- Theoretical and experimental unsteady load distributions 
(first harmonic) on NACA 64A006 airfoil with oscillating flap. 

M = 0.875; k = 0.234. 

OO 


15 


Present calculation 
(M to = 0. 860; a = 0°) 



Figure 4.- Theoretical and experimental steady pressure distributions on 

NACA 64A006 airfoil. 




X 


(a) In-phase component. 



(b) Out -of -phase component. 

Figure 5.- Theoretical and experimental unsteady load distributions 
(first harmonic) on NACA 64A006 airfoil with oscillating flap. 



Case 1 Case 2 Case 3 



(a) Steady pressure. 



(b) Unsteady load distribution. 


Figure 6.- Steady pressure distributions and corresponding unsteady load distri- 
butions on 70-10-13 supercritical airfoil oscillating in pitch at design and 
off -design conditions. 
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Figure 7.- Calculated unsteady load distributions on 70-10-13 supercritical air 
foil oscillating in pitch at off -design conditions. = 0.68; k = 0.10. 





A: kt = 30° 
B: kt = 60° 
C: kt = 90° 


(a) kt = 30°, 60°, 90°. 



(b) kt = 120°, 150°, 180°. 


Figure 8.- Unsteady pressure distribution on 70-10-13 super- 
critical airfoil oscillating in pitch at off -design condition. 
M =0.68; k = 0.10; a=1°+1° sin kt. 

ao 


20 



Nose down I 




G: kt = 210° 
H: kt = 240° 
I: kt = 270° 


! I I I 

-.5 0 .5 1.0 

x 

(c) kt = 210°, 2H0°, 270°. 



kt= 300 C 
kt = 330 C 
kt = 360 ( 


1 I 1 I 

-.5 0 .5 1.0 

X 


(d) kt = 300°, 330°, 360°. 
Figure 8.- Concluded, 



no 

no 



Design (Bauer et a!., ref. 9 ; 

C L = 0. 292) 

Present calculation 

(C L = 0. 308) 


Figure 9.- Steady pressure distribution on 79-03-12 supercritical airfoil in 

design condition. M n 0.79; a = 0°. 

OO 7 









no 

4=r 
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Flat -plate theory 

= 0°), *fa m = 1°) : 79-03-12 airfoil 
= 0°), O' fa m = 1°) : NACA 0012 airfoil 



M 



(a) Unsteady lift coefficient. 

Figure 11.- Magnitude and phase angle (first harmonic) of lift and pitching 
moment on oscillating 79-03-12 supercritical airfoil and NACA 0012 air- 
foil. k = 0.10; a=a m +1° sin kt. 



n 



Flat-plate theory 

• (a = 0°), #(a m = 1°) : 79-03-12 airfoil 
m m 

O (a = 0°), c?(a = 1°) : NACA 0012 airfoil 
m m 



(b) Unsteady pitching -moment coefficient. 


Figure 11.- Concluded. 



20 


79-03-12 Supercritical airfoil 



NACA 0012 airfoil 



(a) M = 0.79; k = 0.10. 

0O 


Figure 12.- Comparison of unsteady load distributions on 79-03-12 supercritical 
airfoil and NACA 0012 airfoil, a = 0° + 1° sin kt. 
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79-03-12 Supercritical airfoil 


a = 0° + 1° sin kt 



AC" 


P. 





NACA 0012 airfoil 



AC" 


P.0 



(b) M = 0.75; k = 0.10. 

GO 


Figure 12.- Continued. 



79-03-12 Supercritical airfoil 



NACA 0012 airfoil 
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(a) kt = 30°, 60°, 90°. 
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F: kt = 180° 


(b) kt = 120°, 150°, 180°. 


Figure 15.- Unsteady pressure distribution on 79-03-12 supercritical 
airfoil oscillating in pitch about midchord axis at = 0.75 
and k = 0.10; a=1°+1° sin kt. 
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(c) kt = 210°, 240°, 270°. 
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(d) kt = 300°, 330°, 360°. 
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Figure 15.- Concluded. 



-.5 


0 

x 


5 


1.0 


(a) b n , n b 1 , 2, 3. (b) a n , n = 1, 2, 3. 

Figure 16.- Chordwise distributions of harmonic components of unsteady pressure 
loadings on 79-03-12 supercritical airfoil oscillating in pitch about mid- 
chord axis at 1^, = 0.75 and k = 0.10; a = 1° + 1° sin kt. 
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Figure 17.- Chordwise distribution of harmonic components of unsteady pressure 
loadings on 79-03-12 supercritical airfoil oscillating in pitch about mid- 
chord axis at M =0.75 and k = 0.30; a = 1° + 1° sin kt. 
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(a) kt = 0°. 


(b) kt = 45°. 


Figure 19.- Unsteady pressure distribution on NACA 0012 oscillating in pitch 
about midchord axis at M =0.79 and k = 0.10; a = 0° + 1° sin kt. 
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Figure 20.- Chordwise distribution of harmonic components of unsteady pressure 
loadings on NACA 0012 airfoil oscillating in pitch about midchord axis at 
M =0.79 and k = 0.10; a=0°+1° sin kt. 
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